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ATTORNEY DOCKET NO. SVG 771 



SYSTEMS AND METHODS FOR CONTROLLING LOCAL ENVIRONMENT 



\ The present invention relates in general to the control of local environment, and 
•y in particular, to systems and methods for controlling the temperature, humidity, 
ij and cleanliness of gas such as air delivered to semiconductor equipment. 

yj{ Background of the Invention 

^ There are many applications in semiconductor manufacturing where tight control 
of temperature, relative humidity, and particles of a process gas such as air 
becomes critical. See Chang & Sze, ULSI Technology (1996), which is hereby 

^ incorporated by reference. 

Foundries, including metrology, lithography, and track areas where photoresist 
and developer is applied to silicon wafers, may require air having a relative 
humidity of 45%+ 0.5% at 24+0.1 C, while other cleanroom areas may require a 
relative humidity of 45%+ 5% at 24+0.5 C. For example, the humidity of air 
delivered to a photoresist spin station must be controlled because it affects the 
evaporation rate of solvent of the resist. To avoid particle contamination, the 
process areas may also require Class 1 (at 0.12 microns) while the surrounding 
areas may only require Class 1000 (at 0.3 microns). The velocity of air is in 
laminar flow region such as 0.35 m/sec. 

Foundries must control other parameters affecting the quality of air, including 
hydrocarbons and other external contaminants, out-gassing of materials 
contacting the air, the pressure of the air, and the electrostatic discharge of the 
system delivering the air to the process. 

Because of these requirements, especially for cleanliness, some foundries seek 
to control the local or mini-environments of critical processes, rather than the 
entire cleanroom, and transport wafers outside mini-environments in SMIF pods 
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1 (standard mechanical interfaces). Chang & Sze, at pages 12-13. The mini- 
environments are cost-effective, but require delivery of precisely controlled gas 
such as air. 

In an air-water system, humidification and dehumidification involves transfer of 
water between a liquid water phase and a fixed air phase which is insoluble in 
the water. Heat and mass transfer effects influence one another so that the 
temperature and relative humidity are coupled together. 

The specific humidity is defined as the mass of water vapor carried by a unit 
mass of dry air. Relative humidity is defined as the ratio of the partial pressure of 
the water vapor to the vapor pressure of the water at the gas temperature, and is 
expressed on a percentage basis so that 100 percent relative humidity means 
saturated air and 0 percent means dry air. 

Accurate regulation of temperature and humidity requires use of closed-loop 
control to compensate for disturbances such as changes in the ambient air 
conditions. The coupling between temperature and relative humidity presents 
challenges to the design of the closed-loop controller. 

The semiconductor industry has attempted to control the temperature and the 
relative humidity of process air in three stages: first, a chiller lowers the 
temperature and humidity of the incoming air below the desired values; second, a 
heater raises the temperature of the air to a desired value, and, third, a steam 
source humidifies the air to achieve the required relative humidity. However, 
boiling water at the steam source is an unstable complex process, and it is 
difficult to provide controlled boiling that generates the required amount of water 
vapor to control humidity within the required tolerances + 0.5% RH. Furthermore, 
the response time of humidity levels resulting from change in the heater power 
supply is too slow and unpredictable. Because the tem perature a nd relativ e 
humidi ty arej soupled and it is difficult to precisely controlling the amount of steam 
used to humidify the air, this conventional approach has significant drawbacks. 
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Summary of the Invention 

The present invention provides a robust and cost-effective alternative to the 
systems used previously because final humidity and temperature can be treated 
as decoupled parameters. 

In a preferred embodiment, a gas-liquid contactor brings liquid into intimate 
contact with gas so the gas leaves the contactor at the saturation temperature, a 
heater heats the gas to a desired temperature, and the absolute humidity of the 
gas is fixed by the inlet liquid temperature to the contactor. 

In a preferred embodiment, the system avoids the difficulty of controlling coupled 
temperature and humidity outputs. In an air-water system, the thermodynamics of 
moist air dictates the dry bulb and saturation temperatures of the air determine its 
relative humidity. With the use of steam tables, the problem of controlling 
temperature and relative humidity can be transformed into one of controlling dry 
bulb and saturation temperature. Because there is no direct coupling in the 
control of these parameters, the invention permits precise control of local 
environments used in the ULSI technology. 
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Brief Description of the Drawings 

Figure 1 is a schematic diagram of one embodiment of the system of the present 
invention. 

Figure 2 illustrates an embodiment of the Direct Contact Cooler-Condenser 
(DCCC) suitable for use in an air-water system. 

Figure 3A is a top view of one embodiment of the sprinkler located in the 
cylindrical housing of the DCCC. 

Figure 3B is a cross-section of a sprinkler arm taken on the line A-A of Figure 3A. 

Figure 3C is a cross-section of a sprinkler arm taken on the line B-B of Figure 3B. 

Figure 3D is a perspective view of a sprinkler arm illustrating the direction and 
spacing of the orifices. 

Figure 3E is cross-section of the sprinkler arm showing its internal structure. 

Figure 4 illustrates an embodiment of the control system suitable for use with the 
present invention. 

Figure 5A-5F are a set of tasks performed by one embodiment of software 
implementing the present invention. 

Figure 6 illustrates an alternative embodiment of the control system. 

Figure 7 is a psychrometric chart showing a process path for humidifying air in 
the system of the present invention. 

Figure 8 is a psychrometric chart showing a process path for dehumidifying air in 
the system of the present invention. 

Figure 9 is a graph illustrating the results of decoupling relative humidity from 
temperature. 
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Detailed Description of the Preferred Embodiments 

I The present invention relates to a system for controlling the local environment of 
t a gas-liquid system. In a preferred embodiment, the system uses a gas-liquid 
^ contactor 100 for humidifying or dehumidifying the gas. As shown in Figure 1 , the 
4 gas-liquid contactor 100 includes a housing filled with a contact media 106, a gas 
/ inlet 102 below the contact media 106, a gas outlet 107 above the contact media 
^ 106, a liquid sprinkler 108 above the contact media 106, and a liquid reservoir 
1 1 1 0 below the housing with an outlet 1 1 2. 

% In a preferred embodiment, the gas-liquid contactor 100 is a direct-contact cooler 
% condenser that operates like a cooling tower, but instead the liquid is heated and 
vo the gas is cooled. The governing equations and solutions are similar. Cooling 
\\ tower theory is dealt with by Mills & Irwin, Heat and Mass Transfer (1995), which 
wis incorporated by reference. 

\^ In a preferred embodiment, the contact media 106 is material with a large surface 
^ area/volume of material, chemically inert to the liquid and gas, which brings the 
^ liquid and gas into intimate contact for sufficient heat and mass transfer to 
\\0 saturate the gas with the liquid after the contact. Some suitable materials are 
y\ tower packing, refrigerated coils, and tubes. McCabe & Smith, Unit Operations 
$(1976 & later editions) and Perry's Handbook of Chemical Engineering (1973 & 
x °vlater editions) provide suitable contact media and are incorporated by reference. 

Conventional cooling tower packing typically has large hydraulic diameters in 
order to reduce pressure drop, but if there is little space, compactness of the 
DCCC 100 may be more important than avoiding greater than typical pressure 
drop. Whatever tower packing is selected, it clearly is preferably to have a large 
surface area per volume. 

The liquid from the outlet 112 enters an inlet 1 16 of a heat exchanger 1 14 
connected to a refrigeration system 146. The chiller, including the refrigeration 
system 146 and the heat exchanger 114, chills the liquid. The liquid returns from 
outlet 1 18 to a pump 120, moves through the flowmeter 140, and enters the gas- 
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liquid contactor 100 at the inlet 104 to the sprinkler 108 forming a closed-loop 
circulation. 

Closed-loop circulation is preferred because the liquid temperature does not 
deviate as much and can be more quickly returned to operating conditions by the 
chiller. Further, closed-loop circulation lowers the expense of purifying water, but 
even so humidification requires make-up and dehumidification removal. Thus, in 
the preferred embodiment, the system provides a make-up inlet 178 and a 
removal outlet 179 at the reservoir 110, and high and low reservoir sensors 174, 
176 indicating when removal or make-up liquid is required. A suitable high and 
low sensor is the NCG Multipoint Level Switch made by Kobold Instruments, Inc. 

A percentage of the liquid in the "closed-loop" is preferably refreshed to meet the 
purity requirements. The capacity of the chiller to bring the liquid to operating 
conditions and the cost of purification sets the appropriate percentage of liquid to 
add to the "closed-loop" through the make-up inlet 178. A preferred percentage 
might be 50% of the circulating liquid per hour. 

It is preferred that the sprinkler 108 distribute liquid uniformly over the contact 
media 106. In one embodiment , the system includes temperature sensors 128 
and 126 to measure the temperature of the liquid entering and leaving the heat 
exchanger 1 14, but the system may in another embodiment include only the gas 
temperature sensors 122 and 144 for control. 

In the embodiment of Figure 1 , a blower 136 moves the gas to a valve 134, a 
positive flow displacement meter (PDFM) 132, and a gas inlet 102 below the 
contact media 106. 

In another embodiment shown in Figure 2, the gas enters inlet 202 from the 
surroundings, the blower 236 is placed downstream of the and valve 

1 34 and flowmeter 1 32 are not used. 

In the preferred embodiment, the gas travels up through the contact media 106 
countercurrent to liquid flowing down. However, it is not strictly necessary to have 
countercurrent flow because any gas-liquid contact achieving the saturated gas 
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at the outlet is satisfactory. After contacting the liquid, the gas leaves at outlet 
107 saturated with vapor and goes to a heater 138 with a heater coil 142 where 
the gas temperature increases to a desired value. A suitable heating coil 142 is a 
Watlow Gordon 1/16 inch diameter cable heater arranged in a snake-like pattern 
in two layers to form a heating grid. As shown in Figure 1 , the system includes 
the gas temperature sensors 124 and 122 to measure the gas entering and 
leaving the gas-liquid contactor 100. 

Thus, the preferred embodiments include two stages. In the first stage, the gas 
passes through contact media where it is cooled and humidified or dehumidified 
by direct contact with the flow of chilled liquid. In this stage, it is important that 
there is sufficient contact media surface area per unit volume to bring the gas 
and the chilled liquid to thermal equilibrium at full saturation by heating the liquid 
and cooling the gas. In the second stage, the gas is heated to the desired dry- 
bulb temperature maintaining its specific humidity unchanged, at which the 
desired relative humidity is achieved because the desired specific humidity was 
achieved in the first stage. From the desired relative humidity and dry-bulb 
temperature, one can extract the corresponding specific humidity and its 
saturation temperature. Saturation temperature can be controlled by controlling 
the temperature of the liquid independently from the dry bulb temperature. 
Similarly, controlling the temperature of the gas in second stage will not impact 
the specific humidity achieved in the first stage. Note the leading digit of the part 
numbers indicates the Figure showing the part in question, while the trailing digits 
indicate like parts are involved. Thus, contact media 106 and 206, shown in 
Figure 1 and Figure 2, respectively, have the same function. 

Figure 2 illustrates a Direct Contact Cooler-Condenser (DCCC) 200 suitable for 
an air-water system. The DCCC 200 preferably includes a cylindrical-shaped 
housing, an air inlet 202, a distributing air space at the bottom of the housing, a 
water inlet 204 and a sprinkler 208 at the top, a reservoir 21 0 with a water outlet 
212 near the bottom, a removal outlet 279 below that, a make-up inlet 278 at the 
top of the reservoir 21 0, and a contact media 206, such as the tower packing 
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mentioned earlier. The housing of the DCCC 200 is preferably chemically inert, 
not off-gassing or exceeding hydrocarbon limits in operation. In one preferred 
embodiment, the internal diameter of the housing of the DCCC 200 is 12 inches 
in diameter and 24 inches in height. The water lines are preferably an inert tough 
material such as !4 inch in diameter Teflon tubing. The inner wall of the reservoir 
210 is larger in diameter than the outer wall of the DCCC 200 forming an inlet 
202 around the circumference of the DCCC 200 to allow uniform airflow into the 
tower. The contact media 206 is tower packing 18 inches in height, leaving an 
open space above of 6 inches. The tower packing is chemically inert, not off- 
gassing, and contains adequate passages to avoid excessive liquid holdup or 
pressure drop. Some suitable tower packings for air-water are tower packing 
media CF650 made by Brentwood Industries, Inc.; black oxidized copper wire 
cloth; cloth of anodized aluminum wire, or Teflon filaments. 

If the contact media 206 is tower packing, it may be held by supporting screens 
(not shown) or held in place by friction of the tower packing against the walls of 
the tower. Preferably, the tower packing has sufficient surface area per unit 
volume to provide good contact between the water and the air for high heat and 
mass transfer rates; and to bring the air and the water into intimate contact 
necessary to achieve thermal equilibrium in a single pass whether parallel or 
countercurrent flow. 

As shown in Figure 2, a non-corrosive chemically compatible pump 220, sitting 
on support 272, withdraws water from the reservoir 210 through the line 212 and 
moves the water through line 216 to a heat exchanger 214. 

One suitable heat exchanger 214 includes a matrix of approximately 100 TECs 
(thermal electric chips) sandwiched between two copper plates such that the cold 
side of the TECs face the same plate and the hot side of the TECs face the 
another plate. Copper tubing with snake-like pattern is brazed on the opposite 
surface of the copper plate that is in contact with the cold surfaces of the TECs. 
Water delivered to the DCCC 200 circulates through the copper tubing. Another 
copper tubing with snake-like pattern is brazed on the opposite surface of the 
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copper plate that is in contact with the hot surfaces of the TECs. Separate water 
runs through this tubing to extract the heat from the TECs. The heat exchanger 
can be remote if desired from the mini-environment to avoid any disturbances. A 
suitable TEC is part no. CP 2.8-32-06L made by Melcor Materials Electronic 
Products Corp. 

The chilled water leaving the heat exchanger 214 at the outlet 218 is delivered 
through the line 204 to the sprinkler 208 which distributes the water on the 
contact media 206. It is important to introduce the water on the contact media 
206 so as not to generate mist, or water particulates. It is also preferred not to 
have any liquid water phase in the air departing the DCCC 200. The air should 
carry water in its vapor-phase only. 

The sprinkler 208 has a set of sprinkler arms 244, 246, and 248, attached to a 
hub 250 capable of free rotation, and preferably, parallel and close to the contact 
media 206 such as one inch or closer, and about four inches or more below a 
gas outlet 207. The relative spacing and the direction of the orifices 399 (Figure 
3) help to prevent entrainment of mist in the air outlet 261 . Thus, the air leaves at 
the outlet 261 saturated and free from water particulates. 

An example of a suitable environment for the embodiment of Figure 2, is an air- 
water system with the following conditions: 

Ambient air temperature: 18-25 C 

Maximum rate of temperature change: + 3 C/hour 

Ambient relative humidity: 35-50% RH 

Maximum rate of relative humidity change: + 5% RH/hour 

Maximum static pressure at output: 1 inch of water 

Temperature control range: 18-24 C + 0.1 C 

Relative humidity: 35 to 50% RH + 0.5% 

Maximum air capacity: 150 CFM 
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Figure 3A is a top view of one embodiment of the rotatable sprinkler 308 in the 
cylindrical tower of the DCCC. As shown, in this embodiment, the sprinkler arms 
344, 346, and 348 attached to a hub 350 are 120 degrees from each other and 
the outer tips of the arms are preferably close, e.g., 1/8 inch from the tower walls 
385. This arrangement when combined with the rotation of the arms promotes 
uniform distribution of the liquid on the tower packing in the cylindrical tower 
(Figures 1-2). It is recognized that other arrangements and number(s) of sprinkler 
arms, such as one or more, e.g., four arms in a cross pattern would be also 
suitable to achieve uniform liquid distribution. Further, it is not necessary that the 
"sprinkler" rotate as it might be a non-rotating shower head, e.g., shaped like a 
disk whose major surfaces are parallel to the contact media. 

Figure 3B is a cross-section of the sprinkler arm 344 taken on line A-A of Figure 
3A. As shown in Figure 3B, an orifice 399 is angled downward, not down or 
sideways, but between, for the following reasons. If the orifice points straight 
down at the contact media it will reduce the mist entrapment in the outlet gas, 
but will exert no component of rotational force on the sprinkler. If, however, the 
orifices 399 are at an angle as shown in Figures 3B, 3D, and 3E, the mist 
entrainment can be prevented, plus the component of reactive force directed 
about the hub can exert rotational force to the sprinkler arm such that the 
sprinkler distributes the liquid uniformly over the contact media without the need 
for a separate motor. 

Figure 3C is a cross-section of a sprinkler arm 344 taken on line B-B of Figure 
3B, showing the tightly spaced orifices 399 which enhance the uniform 
distribution of the liquid on the contact media. 

Figure 3D is a perspective view of the sprinkler arm 344 illustrating the direction 
and spacing of the orifices 399. In a preferred embodiment, the sprinkler 308 has 
a plurality of tightly space orifices 399, e.g., each 1/8 inch apart and 0.050 inches 
in diameter. 
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Figure 3E is a cross-section of the sprinkler arm 344 showing its internal 
structure, including the common passage in the arm 344 for each of the orifices 
399. 

Referring again to Figure 2, the gas leaving at the outlet 261 enters a duct 207. 
The air in the duct 207 travels to a heater 238 close to the exit of the duct 207. 
The heater 238 raises the air temperature to a desired value. The air is pulled 
into a blower 236, leaves duct 207 and enters a diffuser box 265, which has a 
sloped wall 267 such that the air exiting from the diffuser box 267 through a filter 
252 is uniformly distributed on the process below the hood 255. 

Suitable filters for use in an air-water system include a HEPA filter with the 
required efficiency, or an ULPA (Ultra Low Penetration Air) filter having an 
efficiency of 99.9995% or more for particulate diameters greater than 0.12 
microns or another high efficiency air filter. 

The objective of the control system is to regulate the temperature and the relative 
humidity of the gas to user-defined set points, the desired values required by the 
process, and to adjust for disturbances to the system such as changes in the 
temperature of the gas from the cleanroom or another supply source. 

The control systems are now discussed. In the first embodiment, a user sitting a 
computer 464 such as an IBM compatible computer with software described 
below and a conventional analog-to-digital (A/D) converter such as the AT-MIO- 
16DE-10 from National Instruments, which includes software, enters desired 
values for T set point and RH set point, by the keyboard and display monitor. As 
shown in Figure 1, the temperature sensors 122 and 144 measure T saturation 
and T dry bulb of the gas and output analog feedback signals proportional to the 
temperatures to the analog-to-digital (A/D) converter associated with the 
computer 464. 

In one embodiment, the dry bulb temperature and relative humidity have been 
measured after the heater 138 and filter 152 with a combined temperature and 
humidity sensor 144 and 130 such as the HMP233 from Vaisala, Inc. This 
combined sensor 144 and 130 is preferably placed approximately 6 inches 
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\ downstream of the filter 152 outlet. Although a relative humidity feedback signal 
Vis not used in the control, it is entered into the computer 464 for monitoring 
^ purposes. T saturation of the gas out shown in Figure 1 is determined indirectly 
by measuring the temperature of the saturated gas leaving the contactor 100 by 
temperature sensor 122 located just inside the gas-liquid contactor 100. One 
suitable temperature sensor is a resistive-temperature-device (RTD) such as the 
Style RB from Watlow Gordon. 

Figure 4 illustrates the first embodiment of the control system suitable for use 
with the present invention. As shown in Figure 4 t in the first embodiment, the 
control system includes an algorithm consisting of two feedback loops 466, 468, 
and a feedforward element 462 which transforms the T set point and RH set 
point into a corresponding saturation temperature set point 470 sent to 
comparator 460. The comparator 460 takes the difference between set point 470 
and the feedback signal 468, and outputs an error signal to a proportional- 
integral-derivative (PID) compensator 456 which in turn outputs a command 
signal to the chiller 414 to adjust its cooling rate. As is known, a PID takes the 
error signal: (1 ) multiplies it by a first constant; (2) integrates the error signal, and 
multiplies it by a second constant; and (3) takes the error signal's derivative, 
multiplies it by a third constant, and sums (1), (2), and (3), which provides 
accurate and responsive control over a wide range of conditions. See Weber, An 
\ Introduction to Process Dynamics and Control (1973) which discusses PID 
^/controllers in more detail, and is hereby incorporated by reference. 

Similarly, the comparator 458 takes the difference between T set point and 
feedback signal 466 and outputs an error signal to a PID compensator 454, 
which outputs a command signal to the heater 438 to adjust its heating rate. 

^In one embodiment, the transformation to the saturation temperature set point 
470 is achieved by use of empirical 'steam table' equations, which can be found 
in texts such as Irvine & Liley, Steam and Gas Tables with Computer Equations . 
^Academic Press, 1984, which is hereby incorporated by reference. The two 
feedback loops 466, 468, regulate the T dry bulb and T saturation of the outlet 
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gas, while the feedforward 462 ensures T dry bulb and T saturation yield the 
desired relative humidity RH set point defined by the user. 

One embodiment of real-time control software written in the C language is shown 
in the program listing below. The software includes a cooperative multitasking 
scheduler that executes six tasks. In a cooperative multitasking environment, the 
scheduler allocates use of the CPU to each of tasks so each task is given 
processing time and passes back control to scheduler after completion of the 
task. Each task is divided into states which execute in a short time period, so no 
single task monopolizes the CPU, and the tasks appear to be performed 
simultaneously. The software architecture is preferably simple to facilitate 
migration from a PC to an advanced microcontroller such as the Microchip 
PIC17C756 by Microchip Technology, Inc. since all required hardware and 
software can be contained in a compact package. 

Figures 5A-5F show six tasks in state transition logic form performed by one 
embodiment of software implementing the present invention: 

• Task 0 maintains a simple user interface for entering the T and RH set points 
and displaying T dry bulb, RH, and T saturation on the computer monitor and 
permitting shut down of the scheduler if the user presses the escape key on 
the computer. This turns off the chiller 414 and heater 438 and exits the 
control software. 

• Task 1 includes a wait state 507 to initialize the system, including establishing 
communication with the sensors 144 and 130 before initiating control. After 
the system is initialized 509, Task 1 moves into control state 511 where Task 
1 executes the algorithm described in connection with Figure 4, including 
feedforward 462 and feedback loops 466, 468. 

• Task 2 also includes a wait state 51 3 for the same reasons as for wait state 
507. In this embodiment, temperature and humidity sensors 144 and 130 
generate serial data so Task 2 must acquire both serial data and the analog 
signals of sensor 122. The preferred embodiment use only analog sensors 
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because this provides faster sensor response and simplifies the code needed 
in Task 2. 

• Tasks 3 and 4 convert the outputs from heater PID 454 and chiller PID 456 

*\ shown in Figure 4 into pulse-width-modulated (PWM) signals. Horowitz & Hill, 
£ The Art of Electronics (1989) describes PWM signals and is incorporated by 
tg> reference. The PWM signals permit the heater 438 and chiller 414 to interface 
to the computer 464 by conventional solid state relays appropriate for the 
heater 438 and the chiller 414 selected such as those readily available from 
Continental Industries, Inc. or any appropriate power switching device. 
Considerations for selecting solid state relays for the heater and chiller are 
AC versus DC, maximum voltage, maximum current, and switching speed. 
The PWM signals have been generated in software, but a more appropriate 
method would be to use a microcontroller with built-in PWM hardware. 

• Task 5 includes wait state 541 until the user at the computer 464 presses an 
appropriate key to enable the logging of data, T dry bulb, RH, and T 
saturation. Once logging is enabled, Task 5 moves to states 545 where it 
continuously writes the data to the computer hard disk. Another key press 
disables the logging and Task 5 moves back to the wait state. Whether 
control is implemented on a PC or a microcontroller, it is not necessary to 
have Task 5. 

As mentioned previously, the microcontroller embodiment may be preferable to 
replace the personal computer to reduce cost and size. The user interface could 
be reduced to a LCD screen and push buttons. The computational requirements 
of the algorithm have been kept low enough to be performed by inexpensive 
microcontroller such as the Microchip Technology's part no. PIC17C756. 

Figure 6 illustrates a second embodiment of control system suitable for use with 
the present invention. In this embodiment, the system controls T dry bulb and RH 
with two independent feedback loops 666 and 668. A controller 668 and 670 
implement the comparators 658, 660, the compensation of PID 654 and 656, and 
PWM signals driving the heater 638 and chiller 614. Honeywell IAC Division 
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makes suitable single-loop controllers such as the UDC1000 controller. Since 
this control system does not compensate for the coupling of the temperature and 
relative humidity, it is likely to have decreased stability. The advantage of the 
second embodiment is that it can be readily constructed using commercially 
available components. 

The same type of temperature/humidity sensor discussed in the first embodiment 
can be used. The controllers 668 and 670 contain simple functions for setting the 
set points and reading the temperature and humidity sensor. T saturation 
temperature is not used and need not be measured. The outputs of the 
controllers 668 and 670 are PWM signals which interface with the heater 638 and 
the chiller 614 as described in the first embodiment. 

Figure 7, a psychrometric chart, shows humification and relates dry bulb 
temperature of air (abscissa) to specific humidity (ordinate) and % relative 
humidity. As shown in Figure 2 and 7, the air entering inlet 202 of the DCCC 200 
has a relative humidity of 30% at an initial dry bulb temperature of 17 C. This is 
point A shown in Figure 7. As discussed earlier, the air travels through the DCCC 
200, loses sensible heat, and gains vapor from the water contacting it until the air 
is fully saturated. The air leaves the DCCC 200 at outlet 207 at T saturation. This 
is point B shown in Figure 7. The wavy dotted line from A to B indicates that the 
actual temperature-humidity path is not known, nor is it necessary as long as air 
leaves the DCCC 200 at outlet 207 at T saturation, that is, point B. Next, the 
heater 238 heats the outlet air from the DCCC 200 to the desired dry bulb 
temperature T dry bulb and relative humidity RH shown as point C in Figure 7, 
which corresponds to T dry bulb and T saturation shown in the control diagram of 
Figure 4. 

Figure 8, a psychrometric chart, shows dehumidification of air in a system of the 
present invention. As shown in Figure 2 and 8 t the air entering inlet 202 of the 
DCCC 200 has a relative humidity of 50% at a dry bulb temperature of 26 C. This 
is point X shown in Figure 8. As discussed earlier, the inlet air travels through the 
DCCC 200, loses sensible heat, and gains vapor from the water contacting it until 
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the air is fully saturated. The air leaves DCCC 200 at outlet 207 at temperature, T 
saturation. This is point Y shown in Figure 8. The wavy dotted line from X to Y 
indicates that the actual temperature humidity path is not known, nor is it 
necessary as long as the air leaves at outlet 207 at T saturation, that is, point B. 
Next, the heater 238 shown in Figure 1 heats the outlet air from the DCCC 200 to 
the desired dry bulb temperature T dry bulb and relative humidity RH, that is, 
point Z shown in Figure 8. It should be noted that point Z also corresponds to T 
dry bulb and T saturation shown in Figure 4. 

Figure 9 is a graph illustrating the results of decoupling relative humidity from 
temperature of air in the embodiment shown in Figure 2. As shown in Figure 9, 
the user defines RH set point and T set point (see Figures 4 and 6) and 
tolerances of the air at a relative humidity of 45% + 0.5% at 24 +0. 1 C. Although 
T dry bulb, the actual temperature of the air delivered at t = 0 meets the user 
requirements, the actual relative humidity of air, RH, is 46%, and therefore fails to 
meet requirements. After less than about 10 minutes of control, however, the 
system brings relative humidity, RH, of the air within the tolerances established 
by the user. The temperature-humidity path of the air shown in Figure 9 is the 
logged data obtained by execution of Task 5 of the software computer or 
microcontroller shown in Figure 5. 
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The following is a program listing of the software source code in C: 



File: Sched.c 



Bryan Wade 

Sched is a cooperative multi-tasking scheduler-dispatcher. 
Tasks are in the form of state transition logic. Task states are 
entered in the file Tasks.c and scheduling info is entered in this 
file. 



»0){ 



if (task[fgQueue[i]].timeLeft 



#include <stdio.h> 
#include <conio.h> 
#include <time.h> 
#include <dos.h> 
#inctude "sched. h" 
#inctude lasks.h" 
#include "milliclk.h" 

//begin: scheduler ^mmmiimttmmtimmmimmmmmmtm 

#define NUM_TASKS 6 
#define STACK_MAX 2*NUM_TASKS 
struct taskStruct task(NUM_TASKS]={ 
/Tsamp timeLeft nextState priority 
0, 

taskOstateO, 5 



stack(++stacldndex]=fgQueue[i]; 



taskffgQueue[i]].timeLeft 8 taskffgQueueplJ.Tsamp; 
) 

} 

//process the stack 
tf(stacklndexl*EMPTY){ 

//run the task at the top of the stack 

C(task[stack[OD.nextState))0; 

task(stack[0]].timeSinceLast 8 0; 

taskCount[stack[0]]++; 

//update the queue 

for (i=0;i<stackJndex;i++){ 

stackffi=stack[i+1]; 

} 

stacklndex-; 



}else{ 



timeSinceLastV 



{100, 



//run background task 
ff(numBG>0){ 



,0 

{200, 
,0 

{250, 

,0 
{1. 

.0 
0. 

,0 

{2000, 



tasktstateO, 



task2state0, 



task3state0, 
0, 

task4state0, 

0, 
6 



taskSstateO, 



.0 ) 

}; 

//end: scheduler WmitmiMlimmiMlMmMUMmmmmMM 
int StopFlag=0; 
int main(void){ 

int dt; 

int i; 

int stack(STACK_MAXl; 
int stacklndex=EMPTY; 
int fgQueue[NUM_TASKS] ; 
int bgQueue[NUM_TASKS]; 
int bglndex=0; 
int numFG; 
int numBG; 

unsigned long taskCount[NUM_TASKS]; 
unsigned long passes=0; 
unsigned long startTime; 

fillQueues(fgQueue, &numFG, bgQueue, AnumBG, 

taskCount); 

initializationTaskO; 
startTime = 0; 
SetupMilliClockO*. 
while (!StopFlag){ 
//update time stuff 
dt » MitliClock(0); 
passes-*-*-; 

for G=0;i<NUM_TASKS;l++){ 

task[i].tlmeSinceLast +» dt; 

} 

//add states to the stack In order of priority 
if (stacklndex >= STACK_MAX) { 

StopScheduling(STACK_STOP); 

break; 

> 

for(i=0;i<numFG;i++){ 

task(fgQueue[i]].tlmelefl -= dt; 



(*(taskfbgQueue[bg Index]] .nextState)) 0; 

task[bgQueue[bglndexD.timeSinceLast = 0; 

taskCount[bgQueue[bg lndex]]++ ; 

bg Index = 
(bglndex==numBG-1) ? 0 : bglndex+1; 
} 

} 
} 

shutdovvnSchedulerCstartTime.passes.taskCount.stacklndex); 
shutdownTaskO; 

return 0; 
> 

mmmimtmimmummmimmmimmmimmimmmmmii 

void fillQueuesfinl* fgQueue.int* numFG ,int* bgQueue, int* 
numBG, 

unsigned long * taskCount) { 
int i,p; 
int bg=0; 
int fg=0; 

//load the Queues 

for G=0; i<NUM_TASKS; i++){ 

p=task[i]. priority; 

if(p==0){ 

bgQueue[bg] = i; 
bg++; 



}else{ 



} 



fgQueue[p-1] = i; 
fg++; 



> 



taskCount[t]=0; 



*numBG s bg; 
*numFG=fg; 

return; 
) 

IllllilllillliilllllllllllllllilillillllllilllillllillfillillllllllllllliillilU 

void shutdown Scheduler (unsigned long startTime, unsigned 

long passes, 

unsigned long* taskCount, int 

stacklndex){ 

float elapsedTime; 
int i; 

char buffer[20]; 
startTlme^O; 

elapsedTime=(float)((clock0-6tartTime)/CLK_TCK); 
switch (StopFlag){ 

case USER_STOP: 
sprintf (buffer ."User Requested"); break; 
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case STACK.STOP: 
sprintf(buffer,"Task Stack Full"); break; 

case FILE_STOP: 
sprintf(buffer,"File Error"); break; 

case COMM_STOP: 
sprintf(buffer,"Communication Error");break; 
} 

clrscrO; 

prinrffUVt Scheduler Shutdown: %s \n", buffer); 

prinrff\tVt ************************************* **\n"); 

if (ShutdownMilliClockO) 

printf("\t\tTimer ISR terminated 

normally An\n"); 

printfOtVt Scheduling Statistics \n"); 
printfOt\t***************************************Vi"); 

printf("\tMEIapsed Time: 
%8.3f\n",elapsedTime); 

printfOtMScheduler Avg Scan Time: 
%8.3nn",elapsedTime*1000/((float)passes+EPS)); 

for (i=0;i<NUM_TASKS;i++){ 

printf("\t\tTask%d Avg Sample Time: 

%8.3f\n",i, 

elapsedTime*1000/((float)(taskCount[iD+EPS)); 
} 

printfCWTasks Currently in Stack: 
%d\n",stacklndex+1); 
} 

mmtmmmmmimmmmnmmmmmmmmimmmmmm 

void StopScheduling(int stopType){ 
StopFlag=stopType; 

} 

File: Tasks. c 

Bryan Wade 

Tasks is the file where tasks are entered for Sched. Tasks are 
in the form 

of state transition logic. Each state is a function which contains 
ENTRY, 

ACTION and EXIT sections along with other info for the 

scheduler. Any datawtiich needs to be shared between tasks is 

defined as a global variable. 

TaskO - User Interface 

Taskl - Control Algorithm 

Task2 - Read Sensors 

Task3 -Air Heater PWM 

Task4 - Water Chiller PWM 

TaskS - Data Collection 

#include <stdio.h> 
#include <dos.h> 
^include <math.h> 
#include "sched. h" 
#in elude "tasks. h" 
#in elude Iptjo.h" 
#include "screen, h" 
#include "thc_fcns.h" 

/The following 2 files contain code the National Instruments 

data aquistion board. The files are obtained from National 

Instruments.*/ 

#include "nidaq.h" 

#include "nidaqcns.h" 

ffdefine CR OxOO 

//begin: taski global vntimmmimmmmmimmmmmtimmmt 

ffdefine KPTAIR .32 
#defineKITAIR .0015 
ffdefine KDTAIR 0.01 
ffdefine TAIRJLIM .3 
ffdefine KFHEAT1 0.048 
#define KFHEAT2 -0.341 
ffdefine KFHEAT3 6.710 
float TairSetpoint=24; 
float TairError=0; 
float TairErrorlnt=0; 
float TairZeroInt =0; 
float TalrErrorDeriv=0; 




float TairErrorPrev=0; 
float DutyHeaterPredictedsO; 
float DutyHeater=0; 
ffdefine KPTSAT .45 
ffdefine KITSAT .002 
ffdefine KDTSAT 0.01 
ffdefine TSATJLIM .3 
ffdefine KFCHILL1 -0.1156 
ffdefine KFCHILL2 1.5200 
float RHSetpoint=45; 
float TsatSetpoint=0; 
float TsatError=0; 
float TsatErrorlnt=0; 
float TsatZeroInt =0; 
float TsatErrorOeriv=0; 
float TsatErrorPrev=0; 
float DutyChillerPredicted=0; 
float DutyChillen=0; 

//begin: task2 global vaKJIIililllfllliillltlllllilllltlllllUHillllllllltltli 

ffdefine Al DEVICE 1 

ffdefine AfCHAN 0 

ffdefine AI.GAIN 1 

ffdefine TSAT_SLOPE 7.888 

ffdefine TSATJNT -7.500 

Irrt value; 

double vofts; 

float TairMeasured=0; 

float RHMeasured=0; 

float TsatMeasured=0; 

int CommlnffiaIized=FALSE; 

//begin: task3 global wiKJIlilUlilliiliillilllllllilllllUllililliHtlllllliil 

int Outpu1Heater=0; 

float EnableHeater=0; 

float TimerHeater=0; 

float TimeOffHeater=0; 

float PeriodHeater=400; 

//begin: task4 global wKJIIIillillitllilltlliillUillllllliilillllllllllllllll 

int OutputChiller=0; 

float EnableChiller^O; 

float TimerCNIIer=0; 

float TimeOffChiller=0; 

float PeriodChiller=400; 

//begin: tasks global vaisiiilillliltlllllllllllilllllillillfiliillllllliiillill 

ffdefine MAX DATA COUNT 2000 

float CollectData=0;~ 

int DataCount=0; 

float DataTime=0; 

float Data 1 [MAX DATA_COUNT]; 

float Data2[MAX DATA COUNT]; 

float Data3(MAX DATA_COUNTJ; 

float Data4(MAX DATA COUNT]; 

float Data5[MAX DATA COUNT]; 

float Data6[MAX_DATA_COUNTJ; 

float Data7[MAX_DATA_COUNTJ; 

//begin: tasko global s^immmmmmmmmummmmmmmu 

int escKeyHit=0; 
ffdefine NUM_MSG 9 
struct ScreenMsg msg[NUM_MSG]={ 



{ 


"Environmental Control Unit", 


22. 


1}. 


{ 


"OPERATING SETPOINTS", 


4, 


3}, 


{ 


"CONTROL PARAMETERS", 


4, 


7}. 


{ 


"OUTPUT PARAMETERS", 


4. 


12), 


{ 


"DATA COLLECTION", 


4. 


19}. 


{ 


"SENSOR READINGS". 47. 


3}. 




{ 


"INTERMEDIATE VALUES", 


47, 


8}. 


{ 


"OUTPUT VALUES", 47, 


16}. 




{ 


"PRESS [Esq TO EXIT", 


28, 


25} 



}: 

ffdefine NUMJN 9 
struct ScreenTnputttem lnltem[NUMJN]°{ 



TEMP Setpt", 4, 


.05, 


0, 


100, &TairSetpoint}, 






"RH Setpt", 5, 


.25, 


0, 


100, &RHSetpoint}, 






"Zero Tdry", 8, 


1, 


0, 


100, &TairZerolnt}, 






"Zero Tsat", 9, 


1, 


0, 


100, ATsatZeroInt). 






•Enable H",1 3, 1, 


0, 


1. 


AEnableHeater}, 
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"Enable C",14, 


1, 


0. 


1, 


&EnableChiller}, 








-Period H", 16, 


10, 


10, 


5000, 


APeriodHeater), 








"Period C", 17, 


10, 


10, 


5000, 


&PeriodChiller}, 








"Enable", 


20, 


1. 


0, 



1, &CollectData} 

}; 

#defrne NUM_OUT 11 
struct ScreenOutputltem outltem[NUM_OUT]={ 



{ "RHMeas", 5.NULL, &RHMeasured}, 

{ Tair Meas".4,NULL. ATairMeasured}, 

{ Tsat Meas",6.NULL. ATsatMeasured). 

{ Tairlntgr",9,NULL i ATairErrorint}, 

{ "H Duty 1 0.NULL, ADutyHeaterPredicted}, 

{ Tsat$etpt",12.NULL. ATsatSetpoint}, 

{ Tsatlntgr", 13.NULL, ATsatErrorlnt}, 

{ "C Duty * 14.NULL, ADutyChillerPredicted}, 

{ "HDuty",17,NULL. ADutyHeater}, 

{ "C Duty",18,NULL. ADutyChiller}, 

{ "Data Count',20, ADataCount, 
NULL} 



}; 

//begin: taskO s\z\o^liiiiiltlilStililliillllliiiiillilIiiiliillltlIlllHlllill 

void task0state00{ static cont = FALSE; 
if (cont==FALSE){ 

cont=TRUE; 

//begin: ENTRY function 

SetupScreen(msg,NUM_MSG, 

inltem,NUM_IN, outltem.NUM_OUT, 

AescKeyHit, NULL); 

//end: ENTRY function 

> 

//begin: ACTION function 
DoScreenlnputO; 
DoScreenOutputO ; 
//end: ACTION function 
if(escKeyHit){ 

//begin: EXIT function 

StopScheduling(USER_STOP); 

//end: EXIT function 

cont=FALSE; 

ta$k[0].nextState « taskOstateO; 

} 

if (!Commlnitiali2ed){ 

//begin: EXIT function 
StopScheduling(COMM_STOP); 
//end: EXIT FUNCTION 
cont= FALSE; 

task[0].nextState « taskOstateO; 

) 

} 

//begin: taski ^x^mmimmimmmmmmmmmmmmnmii 

void taski stateOOC static cont * FALSE; 

lf(cont==FALSE){ 
cont=TRUE; 
//begin: ENTRY function 
//end: ENTRY function 

} 

//begin: ACTION function 
//end: ACTION function 
if (Commlnitialized){ 

//begin: EXIT function 

//end: EXIT function 

cont=FALSE; 

task|1].nextState = tasklstatel; 

} 

} 

void tasklstatel 0( static cont ° FALSE; 

if (cont==FALSE){ 
cont°TRUE; 
//begin: ENTRY function 
//end: ENTRY function 

} 

//begin: ACTION function 

/// DRY BULB TEMPERATURE CONTROL: PID w/ 

FF// 

/•Compute P term */ 




TairError n TairS etpoi nt-Ta iiMeasu red ; 

/•Compute I term*/ 

TairErrortnt +* 
TalrError*(noa1)(task(1].timeSinceLast)/1000; 

if (TairZerolnt || fabs(TairError)>TAIRJLIM){ 
TairErrortnM).; 

} 

rCompute D term*/ 
TairErrorDeriv = (TairError - 
TairErTorPrev)/((float)(task{1].timeSinceLast)/1000); 
TairErrorPrev = TairError; 

rCompute Feedforward*/ 
DutyHeaterPredicted = KFHEAT1*(TairSetpoint- 
(1 .+KFHEAT2)*TsatMeasured-KFHEAT3); 

/•Compute Output*/ 
DutyHeater = KPTAIR*TairError + 
KITAIRTairErrorlnt ♦ KDTAIRTairErrorDeriv 

+ DutyHeaterPredicted; 
if (DutyHeater<0) DutyHeaten=0.; 
else if (DutyHeater>1.) DutyHeater=1.; 

//////// SAT TEMPERATURE CONTROL: PID w/ FF 

/// 

/•Determine 

Tsat Setpoint */ 

TsatSetpoint = 

Tsat_steamTable(RHSetpoint/100*Psat_steamTable(TairSetpo 
int+273.16))-273.16; 

rCompute P 

term*/ 

TsatError = TsatSetpoint - TsatMeasured; 

^Compute I term 

7 

TsatErrorlnt += 
TsatError*(float)(task(1].timeSinceLast)/1 000; 

if (TsatZeroInt || fabs(rsatError)>TSAT_ILIM){ 
TsatErrorinM).; 

} 

/"Compute D 

term*/ 

TsatErrorDeriv = (TsatError - 
TsatErrorPrev)/((float)(task[1].timeSinceLast)/1000); 
TsatErrorPrev « TsatError; 

/"Compute 

Feedforward */ 

DutyChillerPredicted « KFCHILL1*TsatSetpoint + 

KFCHILL2; 

/"Compute 

Output */ 

DutyChiller ° KPTSATTsat Error + 
KITSATTsatErrorlnt + KDTSATTsatErrorDeriv 

+ DutyChillerPredicted; 
if (DutyChillerO) DutyChiller=0.; 
else if (DutyChiller>1.) DutyChiller=1.0; 
//end: ACTION function 
if (FALSE) 
{ 

//begin: EXIT function 
//end: EXIT function 
cont=FALSE; 

task[1].nextState ° tasklstatel; 

} 

) 

//begin: task2 %\^\^mmtmmmmmimmumnmmmmmmm 

void task2state0(){ static cont ° FALSE; 
if(cont==FALSE){ 

cont=TRUE; 

//begin: ENTRY function 

//end: ENTRY function 

) 

//begin: ACTION function 
//end: ACTION function 
rf (Commlnitialized){ 

//begin: EXIT function 

//end: EXIT function 

cont=FALSE; 

task(2].nextState ° task2state1; 

> 

} 

void task2state1 Q{ static cont ° FALSE; 
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if (cont==FALSE){ 

cont=TRUE; 

//begin: ENTRY function 

txCharCS*); 

txCharCEO; 

txCharfN*); 

txCharfDO; 

txChar(CR); 

//end: ENTRY function 

} 

//begin: ACTION function 
//end: ACTION function 
if (TRUE){ 

//begin: EXIT function 
if 

(getHMP230Readings(&RHMeasured t &TairMeasured» { 

StopScheduling(COMM_STOP); 
} 

//end: EXIT function 
cont= FALSE; 

task[2]. next State = task2state2; 

} 

} 

void task2state20{ static cont = FALSE; 
if(cont==FALSE){ 

cont=TRUE; 

//begin: ENTRY function 

AI_Clear(AI DEVICE); 



ue,& volts); 
TSATJNT; 



ALRead(At_DEVICEALCHAN,AI_GAIN 1 &value); 
ALVScale(AI_DEVICE,ALCHAN^t_GAiN,1,0.0,val 
TsatMeasured = TSATSLOPE*volts + 
//end: ENTRY function 



) 

//begin: ACTION function 
//end: ACTION function 
if (TRUE){ 

//begin: EXIT function 
//end: EXIT function 
cont=FALSE; 

task[2].nextState = task2state1; 

> 

} 

//begin: task3 s\atestMmmmmmmtMW!WWWW! 
void task3stateO(){ static cont = FALSE; 
if(cont==FALSE){ 

cont=TRUE; 

//begin: ENTRY function 

Pin2lo(); 

OutputHeater=0; 
//end: ENTRY function 

} 

//begin: ACTION function 

//end: ACTION function 

if (EnableHeater==1.){ 

//begin: EXIT function 
//end: EXIT function 
cont=FALSE; 

task[3].nextState = task3state1; 

} 

} 

void task3state10{ static cont ° FALSE; 
if(cont«=FALSE){ 

cont=TRUE; 

//begin: ENTRY function 

TlmerHeater=0; 

TimeOfTHeater=DutyHeater*PeriodHeater; 
Pin2hl0; 

OutputHeater°1; 
//end: ENTRY function 

} 

//begin: ACTION function 

TimerHeater+a (float)task[3] .timeSinceLast; 

//end: ACTION function 



if (TimerHeater>TirneOffHeater && 
DutyHeater>=1.){ 

//begin: EXIT function 
//end: EXIT function 
cont*FALSE; 

task[3].nextState * task3state1; 

) 

if CTimertteater>TimeOffHeater && DutyHeater<1 .){ 
//begin: EXIT function 
//end: EXIT function 
cont=FALSE; 

taskpj.nextState « task3state2; 

} 

if (EnableHeater<1){ 

//begin: EXIT function 
//end: EXIT function 
cont=FALSE; 

taskpl.nextState = task3state0; 

} 

} 

void task3state2(){ static cont = FALSE; 
if(cont=FALSE){ 

cont=TRUE; 

//begin: ENTRY function 

Pin2lo0; 

OutputHeater=0; 

//end: ENTRY function 

} 

//begin: ACTION function 

TimerHeater+= (float)taskf3].timeSinceLast; 

//end: ACTION function 

if (TimefHeater>PeriodHeater && DutyH eater <=0){ 
//begin: EXIT function 
//end: EXIT function 
cont= FALSE; 

task[3].nextState= task3state2; 

} 

if (TimerHeater>PeriodHeater && DutyHeater>0){ 
//begin: EXIT function 
//end: EXIT function 
cont=FALSE; 

task[3]. next State = task3state1; 

) 

> 

//begin: task4 ^X^HIHUmHUIUIIflUIIIIHIUniUilllUniHIIIIlin 
void task4state00( static cont = FALSE; 
if (cont==FALSE){ 

cont=TRUE; 

//begin: ENTRY function 

Pin3lo0; 

OutputChiller=0; 

//end: ENTRY function 

) 

//begin: ACTION function 
//end: ACTION function 
if (EnableChiller==1.){ 

//begin: EXIT function 

//end: EXIT function 

cont= FALSE; 

task[4).nextState = task4state1; 

} 

} 

void task4state10{ static cont ° FALSE; 
if(cont==FALSE){ 

cont=TRUE; 

//begin: ENTRY function 

TimerChiller=0; 

TlmeOffChilIep=DutyChiller*PeriodChiller; 
Pin3hi0; 
OutputChiller=1; 
//end: ENTRY function 

) 

//begin: ACTION function 

TimerChiUer+= (fl oat)t as <4] Jim ©Since Last; 

//end: ACTION function 

if CTimerChiller>TimeOffChiller && DutyChiller>=1 .){ 
//begin: EXIT function 
//end: EXIT function 
cont=FALSE; 
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task(4].nextState = task4state1 ; 

} 

if (TimerChiller>TtmeOffChil]er && DutyChiller<1.){ 
//begin: EXIT function 
//end; EXIT function 
cont= FALSE; 

task(4].nextState = task4state2; 

} 

if(EnableChitler<1){ 

//begin: EXIT function 
//end: EXIT function 
cont=FALSE; 

task[4].nextState = task4state0; 

} 

} 

void task4state20{ static cont = FALSE; 
if(cont==FALSE){ 

cont=TRUE; 

//begin: ENTRY function 

Pin3lo(); 

OutputChiller=0; 

//end: ENTRY function 

} 

//begin: ACTION function 

T1merChiiler+B (float)task[4].timeSinceLast; 

//end: ACTION function 

if (TimerChiiler>PeriodChil]er && DutvChiller<=0){ 
//begin: EXIT function 
//end: EXIT function 
cont=FALSE; 

task[4].nextState * task4state2; 

) 

if (TimerChiller>PeriodChiller && DutyChiiter>0){ 
//begin: EXIT function 
//end: EXIT function 
cont= FALSE; 

task[4].nextState = task4state1; 

} 

} 

//begin: taskS $\aiesJilillllliil!l!ili!UIlllllllllllHIIIiIII!lii!llilIllillil 
void task5state00{ static cont = FALSE; 
if (cont==FALSE){ 

cont=TRUE; 

//begin: ENTRY function 

DataCount=0; 

//end: ENTRY function 

} 

//begin: ACTION function 

DataCount=0; 

//end: ACTION function 

if (CollectData==TRUE){ 

//begin: EXIT function 

//end: EXIT function 

cont=FALSE; 

task(5].nextState = task5state1 ; 

} 

} 

void task5state10{ static cont » FALSE; 
if (cont==FALSE){ 

cont*TRUE; 

//begin: ENTRY function 

DataTime=0; 

//end: ENTRY function 

} 

//begin: ACTION function 

DataTime+= (f1oat)(task(5].timeSinceLast)/1 000.; 

if (DataCount<MAX_DATA_COUNT) { 

Data1[DataCount]»DataTlme; 

Data2[DataCount]=TairMeasured; 

Data3[DataCount]=TairSet point; 

Data4(DataCountj»TsatMeasured ; 

Data5[DataCount]=TsatSetpolnt; 

Data6[DataCount]=RHMeasured; 

Data7[DataCount++]=RHSetpoint; 

) 

//end: ACTION function 
if (CollectData==FALSE){ 

//begin: EXIT function 

//end: EXIT function 

cont=FALSE; 



} 



task[5].nextState * taskSstateO; 



} 

//begin: \nma\aa\\onTaBk/mmfmw/m/wmm//W/m/W/ 
void initlalizationTaskO { 

//initialize serial comm 

inKCommunicationsO; 

Commlnitiali2ed=TRUE; 

//initialize DAQ 

USE_E_Series_AIO; 

AI_Clear(AI_DEVICE); 

//initialize Screen 

clrscrO; 

) 

//begin: ^m^i^^wmimmmmmimmmimmmmmmi 

void shutdownTaskO { 
infix; 
FILE *fp; 

//Set all outputs low 

Pin2lo0; 

Pin3lo0; 

//Write data to file 

fp=fopenfdata.txt",\vO; 

for (ix=0; ix<DataCount; ix++) { 

fprintf(fpr%fVt%f\!%f\t%nt%f\t%f\t%f\n" l Data1 px],Dat 
a2[ix],Data30x], 

Data4[ix] 1 Data5[ix] J Data6[ix],Data7[ixD; 
} 

fclose(fp); 

} 



File: Screen.c 

Bryan Wade 

Creates a simple text ui for real-time systems. The interface 
entities are 1) static messages (e.g. a title), 2) outputs with 
labels, 

3) inputs with labels (increment/decrement by user-defined 
amount), and 4) 
access to the ESC key. 
Use: 

1) #include "screen.h" 

2) define struct arrays: ScreenMsgQ .Screen ItemQ (inputs) 
ScreenttemQ (outputs) 

3) void SetupScreen(struct ScreenMsg* msg, int numMsgs, 
struct Screenlnputltem input, 

int numtnputs, struct ScreenOutputttem* 
output, int numOutpus, int* ESC_flag); 

4) void DoScreenlnputsO; - 
process keyboard buffer for input commands 

5) void DoScreenOutputsO; - 
updates outputs (one output per call) 

#lnclude <conio.h> 

ftinclude <stdio.h> 

#include "screen.h" 

static struct ScreenMsg* msg; 

static int maxMsglndex; 

static struct Screenlnputltem* in item; 

static int maxlnputlndex; 

static struct ScreenOutputttem* outltem; 

static int maxOutputlndex; 

static int* escFlag; 

static int* spaceFlag; 

mmmmmmnmmmmmimmmmnmmm 

void DoScreenlnputO{ 
intc; 

static int index°0; 
while(kbhitO){ 

c » getchO; 

switch (c){ 

case UP: 

gotoxy(INPUT_CURS_COL,inKempndex].pos); 
printff "); 
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index = (index ==0)? maxlnputlndex : index-1; 

gotoxy(INPUT_CURS_COL,inltempndex].pos); 

printff>"); 

break; 

case DOWN: 
gotoxy(INPUT_CURS_COL.inHemPndex].pos); 
printff "); 

index = (index ==maxlnputlndex)? 0 : index+1; 

gotoxy(INPUT_CURS_COLJnKem[index].pos); 

printff>*); 

break; 

case LEFT: 

*(inttempndex].address)-= inltempndexj.defta; 
if (*(inltempndex].address) < inttempndex].min){ 

*(inltempndex].address) n 

inltempndex].min;} 

gotoxy(INPUT_VALU_COL,inltemPndex].pos); 
printff%7.2f V(inltempndex].address)); 
break; 

RIGHT: 



*(inltemfindox].address)+=inttempndex].detta; 
if (*(inltempndex].address) > inltempndex].max){ 

*(inttem[index]. address) = 

inKempndex].nriax;} 

gotoxy(INPUT„VALU_COL,inltempndex].pos); 

printff %7.2f ",*(inltempndex].address)); 

break; 



»1; 



break; 



case ESC: 



case SPACE: 



*spaceFlag=1; 
break; 

} 

} 

} 

miummimiimmmiMimmmmmmtiummutmiMi 

void DoScreenOutput(){ 

static int index=0; 

gotoxy(OUTPUT_VALU_COL,outltempndex].pos); 
if (outltempndex].iPtr!=NULL) 

printff%d V(outltempndex].iPtr)); 

else 

printff%7.2f V(outltempndexJ.fPtr)); 
index = (index<maxOutput Index) 7 index +1:0; 

} 

mmmtmmimmmmmmmmmmmimmmmmnmm 

void SetupScreen(strud ScreenMsg* m, Int nm, struct 
Screen In put Item* li.fnt ni, 

struct ScreenOutputltem* oi.lnt no, int* e, Int* spc){ 
int i; 

//set the global ptrs 
msg=m; 

maxMsg Index = nm-1; 
inltem=H; 

maxlnputlndex=ni-1; 



outHem=oi; 

maxOutputlndex=no-1 ; 

escFlag = e; 

spaceFlag * spc; 

//setup messages 

for (i«0; l<*maxMsglndex; !++){ 

gotoxy(msgP]Jt.msgp].y); 

printff %s" t msg[i].text); 

} 

//setup inputs 

for (i=0; i<=maxlnputlndex; i++){ 

gotoxy(INPUT_r4AME_COL,inltem[i].pos); 
printff %s",inltemp].name); 

gotoxy(INPUT_VALU_COL,inltemp].pos); 

printff %7.2T ) *(inltem[i] .address)); 

> 

gotoxy(INPUT_CURS_COL,inltem[0].pos); 

printff>-); 

//setup outputs 

for (i=0; i<=maxOutpu1lndex; H-+){ 

gotoxy(OUTPUT_NAME_COL,outHemp].pos); 
printff %s",outltemp].name) ; 

gotoxy(OUTPUT VALU COL,outttem[i].pos); 
if(outitemplJPtr!»=NULL) 
printff%d **(outltemp].iPtr)); 

else 

printff%7.2f Y(outltemp].fPtr)); 



File: The fensx; 



Bryan Wade 

This file contains functions for the DCCC controller. 
This includes the steam table equations and routines for 
communicating with 
the HMP230 sensor. 



"7 



tfinclude "thc_fcns.h" 

#include <math.h> 

#include "bios.h" 

#include <string.h> 

#defineCOM0 //0=COM1, 1=COM2 

#define TRUE 1 

#define FALSE 0 

#define CR OxOD 

#define NULO 

#define READY TO WRITE 0x2000 
#define DATA_READY 0x100 

/The following equations are based on material in the text: 
Steam and Gas Tables with Computer Control_, Irvine and 
Liley, 1984. 
7 

float 

Tsat_steamTable(float P){ 

const double A = ,426776e2; 
const double B » -.389270e4; 
const double C a -.9486 5e1; 



return A + B/(log(P) + C); 
} 

float 

Psat_steamTable(float T){ 
const float AO a , 
const float A1 = • 
const float A2 * ■ 
const float A3 ° . 
const float A4 a - 
const float AS = 
const float A6 = . 
const float A7 ° • 
const float A8 ° . 
const float A9 D . 
const float A10=» 



,104592e2; 

..404897e-2; 

..417520e-4; 

.36851 0e-6; 

..101520e-8; 

.86531 0e-1 2; 

903668e-15; 

.1996900-17; 

.779287e-21; 

.191482e-24; 

•.396806e4; 
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const float A11° .395735e2; 

return 

exp(AO+A1T-+A2*D(Jw(T ( 20+A3*powa.3.)+A4 # powCT,4.)+A5*p 
ow(T,5.)+ 

A6*pow(T,6.)+A7*pow(T,7.) + A8*pow(T.8.)+A9*povv( 
T,9.)+A10/(T-A11)); 

} 

void 

initCommunications(void){ 

int init_data, doneFlag; 

charch, bufI10]; 

unsigned int ix=0,tries; 

init_data ° (_COM CHR8 |_COM STOP1 | 
_COM_NOPARITY |_COM_9600); 

_bios_serialcom(_COMJNIT, COM, inrt_data ); 

//Repeat request until get a good response 

//(Some times HMP230 takes a couple tries to work 

at first) 

for (doneFlag=FALSE;doneFlag==TRUE;) { 
txChar(CR); 
txCharCS*); 
txCharCE'); 
txCharCN*); 
txCharCD'); 
txChar(CR); 

for (ix=0,tries=0; ix<11;) { 
ch=rxChar(); 

tf(ch!=NUL)bufpx++]=ch; 
if (++tries>20000) break; 

} 

if (buf[3K»7 && bufl3]==\') 

doneFlag»TRUE; 
} 

} 

char 

rxChar(void){ 

char ch; 

unsigned int status; 

status = _bios_serialcom(_COM_STATUS, COM, 

0); 

if (status & DATA READY){ 
ch = 0x7f & 

J>ios_seriatomC_COM_RECEIVE. COM, 0); 

retum(ch); 

} 

retum(NUL); 
} 

void 

txChar(char ch)( 

unsigned int status; 
/Await 

status = READY TO_WRITE & _bios_serialcom( 
COM STATUS, COM, 0); 

while ( status != READY_TO WRfTE && 0){ 
status = READY TO WRITE & 
_bios_serialcom CCOM_STATUS. COM, 0); 
) 

//send char 

_bios serialcom( COM SEND, COM, ch ); 

} 

int 

getHMP230Readings(f1oat *rh, float *t) { 
char ch, buf[10]; 
unsigned int ix=0,tries; 

for (lx=0,tries=0; ix<6;) { 

ch»rxChar(); 

rf(chl°NUL){ 

buf[ix++]«ch; 

} 

rf(++tries>40000){ 

txChar(CR); 
return 0; 

} 

} 

bufpx)=NUL; 

if (buf|3) '.•) *rh* atof(buf); 
else 

return 0; 
for (bc=0,tries=0; ix<5;) { 
ch=rxCharQ; 



if(ch!=NUL){ 

bufpx++]=ch; 

} 

if(++tries>40000){ 

txChar(CR); 
return 0; 

) 

) 

buf[ix}=NUL; 

if(buf[2]™7) *t= atof(buf); 

return 0; 

^ # 

File: Milltelk.c 



Bryan Wade 
Milliclk uses the timer interrupt to maintain alms 
resolution clock. 
Use: 

1) ffinclude "millidk.h" 

2) void SetupMilliClockO; sets up timer interrupt and resets 
clock to 0 

3) int MilliClock(l); gets time (ms) since last reset 

4) int MilliClock(O); gets time (ms) since last reset & resets 
clock to 0 

5) void ShutdownMitliClockO; resets interrupt to DOS 
default - very important! 

ffinclude <dos.h> 
include -millidk.h" 
void interrupt (*old_vect)0; 
static int count; 
static int milliCount; 
static int running=0; 

mmimumimmmmmmimmmimmmmmmmimmmm 

void SetupMiiliClockO{ 
if (Irunning) { 

milliCount = 0; 

old_vect=getvect(CK_VECT); 
disableO; 

setvect(CKJ/ECT, timerJSR); 
setalarm(1 .0); 
enabteO; 
running 3 1 ; 

) 

} 

nmmmimmmmmimttmmmmmtttimmmnmmmmim 

intShutdownMilliClockO{ 
If (running){ 

disableO; 

setalarm(0.0)/*set back to DOS 

default*/ 

setvect(CK_VECT,old_vect); 

enableO; 

running=0; 

} 

retum(1); 
) 

llllllllllliilitiitlllllliliHitlililliiillllHiilllllllilllllHIillllltililiili 
int MilliClock(int flag){ 

int temp * milliCount;; 

if (flag==0) milliCount ° 0; 
return temp; 
} 

immmiimumtmimimmmimmmmmiiimmmmmmm 

void interrupt timer_ISR(void){ 

/"manipulate floating point processor stack*/ 

asm sub SP.94 

asm mov BP.SP 

asm fsave(BP] 

asm fwait 

mHUCount++; 

Tsend EOIV 

outport(0x20,0x60); 

/•manipulate floating point processor stack*/ 
asm mov BP.SP 
asm frstorfBP] 
asm add SP.94 
asm fwaft 
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} 

mmmtimtmmHmmmmimmmmimmmmmmmmmm 

int setalarm (double ms){ 

if((ms>0.0) && (ms<=MAX_MS)){ 

count=(CLOCKFREQ*ms) / 1000.0 + 

0.5; 

outportb(TIMER__CTL, 0x36); 

outportb(TIMER, count&OxFF); 
outportb(TIMER, (count»8)&0xFF); 
}else if (ms<=0.0){ 

outportb(TIMER CTL, 0x36); 
outportbCTIMER. DOS_COUNT&0xFF); 
outportb(TIMER. 
(DOS COUNT»8)&0xFF); 

count=DOS COUNT; 

}else{ 

retum(1);/*error return*/ 

} 

retum(0)/*normal return*/ 
} xxi ^ x 

File: LPTJO.c Bryan Wade 

Functions for using parallel port for digital I/O. Now that we 
have the N I data acquisition board it would be easier to use it 
for I/O rather than this. 



#include <dos.h> 

tfdefine BASE 0x378 /* Usually Ipt1 is 0x378 & Ipt2 is 0x278. */ 
static int dataByte = 0x00; 
static int inputByte » 0x00; 
static int controlByte « 0x00; 
void Pin2hi0{ 

dataByte |=0x01; 

outportb(BASE+0,dataByte); 

) 

void Pin2lo0{ 

dataByte &= OxFE; 
outportb(BASE+0,dataByte); 

} 

void Pin3hi0{ 

dataByte |° 0x02; 
outportb(BASE+0,dataByte); 

) 

void Pin3lo0{ 

dataByte &« OxFD; 
outportb(BASE+0,dataByte); 

) 

void Pin4hi0{ 

dataByte |« 0x04; 
outportb(BASE+O t dataByte); 

> 

void Pin41o0{ 

dataByte &= OxFB; 
outportb(BASE+0.dataByte); 

} 

void Pin5lo0{ //inverted 

dataByte |= 0x08; 
outportb(BASE+0,dataByte); 

) 

void Pin5hi0{ //inverted 

dataByte &= 0xF7; 
outportb(BASE+0,dataByte); 

) 

void Pin6hi0{ 

dataByte |= 0x10; 
outportb(BASE+0,dataByte); 

) 

void Pin6lo0{ 

dataByte &° OxEF; 

out port b(B AS E+0, dataByte); 

) 

Int ShutdownLPTO{ 

outportb(BASE+0,0x00); 

retum(1); 
} 

File: Sched.h* 

//begin: sched ^x^mmnmiimmtmimmummmimmmmm 



tfdefine EMPTY -1 
#define EPS 1e-9 
tfdefine USER.STOP 1 
#define STACK_STOP 2 
tfdefine FILE_STOP 3 
#defineSTEPPING_STOP 4 
#define COMM_STOP 5 
struct taskStruct{ 

int Tsamp; 

int timeLett; 

void (*nextState)0; 

int priority; 

int timeSinceLast; 

); 

//begin: sched function ^o\o^esJliilll!HillHiillilllillliiiiliiillillli!fl 
void fillQueuesOnt* fgQueue.int* numFG,int* bgQueue.int* 
numBO, 

unsigned long* taskCount); 
void shutdownScheduler( unsigned long startTime, unsigned 
long passes, 

unsigned long* taskCount, int stacklndex); 
void StopScheduling(int stopType); 



File: Tasks.h 

//taskO states prototypes 

void taskOstateOO; //user interface 

//taskl state prototypes 

void tasklstateOO; //goon 

void taskl state! 0; //process control algorithm 

//task2 state prototypes 

void task2state00; //goon 

void task2state1 0; //request RHJair and wait for response 

void task2state20; //readTsat 

//task3 state prototypes 

void task3state00; //heater output disabled 

voidtask3state10; //heater PWM hi 

void task3state20; //heater PWM lo 

//task4 state prototypes 

void task4state00; //Chiller output disabled 

void task4state1 Q; //Chiller PWM hi 

voidtask4state20; //Chiller PWM lo 

//taskS state prototypes 

void taskSstateOO; //Don't collect data 

voidtask5state10; //Collect data 

//inttalizationTask prototype 

void initial izationTask(void); 

//shutdownTask prototype 

void shutdownTask(void); 

//end: state prototyping///////////////// 

//begin: Tasks module defines 

extern struct taskStruct taskfl; 

tfdefine FALSE 0 

tfdefine TRUE 1 

tfdefine HUGE1.0e9 



File Screen.h 

tfdefine UP '8' 

«define DOWN '2' 

•define LEFT '4' 

•define RIGHT *& 

•define ESC 0x1 b 

tfdefine SPACE" 

tfdefine INPUT NAME COL 6 

tfdefine INPUT_CURS COL 5 

tfdefine INPUT VALU COL 18 

tfdefine OUTPUT NAME COL 49 

tfdefine OUTPUT_VALU_COL 61 

void DoScreenlnput(void); 

void DoScreenOutput(void); 

void SetupScreen(struct ScreenMsg*,int, struct 

ScreenlnputHem*,int, 

struct ScreenOutputltem*,int,int*,int*); 
struct Screenlnputltem { 
char name[12]; 
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Int pos; 
float delta; 
float min; 
float max; 
float* address; 

}; 

struct ScreenOutputltem { 
char name[12]; 
int pos; 



# 



int* IPtr; 
float* fPtr; 

}; 

struct ScreenMsg { 

char text [80] ; 

intx; 

inty; 

}; 

Typedef unsigned char BYTE; 
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